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ABSTRACT 

We report observations of ionized and warm molecular gas in the extended regions of the 
central galaxies in several cooling flow clusters. These show that both gas phases are present 
in these clusters to large radii. We trace both Paa and H2 lines to radii in excess of 20 kpc. The 
surface brightness profiles of the two phases trace each other closely. Apart from very close 
to the central AGN, line ratios in and between the phases vary only slowly with position. 
The kinematics of the phases are indistinguishable and away from the influence of the central 
AGN both the mean and dispersion in velocity are low (< 100km s _1 ) ruling out kinematic 
support of the molecular gas. All of the above indicates that the mechanisms for heating the 
molecular gas and ionizing the HII regions are highly coupled. The highest surface brightness 
emission within a few kpc of the central AGN is distinct, both kinematically and thermally 
from that at larger radii. The relative strengths of the Paa to the H2 lines indicates a source 
of UV excitation rich in EUV relative to FUV photons, e.g. a black body with a temperature 
> 10 5 K. 
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1 INTRODUCTION 



The density of the hot (10 7 - 10 8 A") intracluster medium (ICM) at 
the centres of many galaxy clusters is high enough for it to cool 
on timescales less than the age of the cl uster, potentia lly lead- 
ing to the formation of a cooling flow ("see lFabian [l994 for a re- 
view). Recent observations of galaxy clusters with Chandra and 
XMM-Newton have shown that the situation at the centres of clus- 
ters thought to contain cooling flows is considerably more com- 
plex than previously thought. Chandra and XMM-Newton data in- 
dicates that the size of a cooling region is typically smaller than 
was indicated by earlier ROSAT data fe.e. lPeres et alll99a) . cool- 
ing may be intermittent, and there is considerable interaction be- 
tween the hot phase at the centre of a clust er and any ra dio plasma 
emitted by the central radio so urce iFabian et alIl20oiLlBoh ringer. 
l200ltlMcNamara et all 1200 lh . Chandra and XMM-Newton spec- 
troscopy of the central regions of clusters reveals a lack of line 
emission expected from gas at ~ 2 keV, implying that the amount 
of gas cooling through these temperatures is considerably less than 
the mass deposition rates determined from ROSAT data (e.g. see 
IPeterson et all2003l) . 

The existence of atomic and ionised gas emitting at 10 4 K in 
and around the dominant galaxies in many of th ese cooling flow 
clusters has been kn own for two decades (e.g. iHu et alj. Il983t 
Heck man et all ll998). Despite clearly being statistically correlated 
with the presence of the cooling flow in the cluster, the emission- 
line gas discovered in these early studies has too high a surface 
brightness and luminosity to be simple recombination radiation 
from material cooling through 10 4 K. It has long been recognised 



that this relatively high surface brightness emission must be ex- 
cited by one or more mechanisms: hot stars, shocks, turbulence or 
X-ray heating. 

Cooler molecular material associated with the optical nebu- 
losity has recently been detected from its K-band H2 line emission 
tfaffe & Bremed. 1 19971 : iFalcke et all Il998l : IWilman et all l200(j 
Donah ue et all l2000ti . This material is dense (n e > 10 5 cm -3 ) 
and it s occurance corr elates with that of atomic nebular emis- 
sion IJaffeetallEoOll) . This, plus detection of large amounts of 
cold «50K) CO in cooling flows lEdsej <200ll) indicates that large 
amounts of cool material is present at the centres of these clusters, 
even if the new X-ray data indicates weaker cooling than was pre- 
viously estimated. 

With the exception of the CO observations whose spatial reso- 
lution is so far quite poor, most previous studies of gas in brightest 
cluster galaxies (BCGs) have only detected gas within a few kpc 
of the nuclei, where the energetics of th e gas are w ithout doubt 
dominated by the AGN activity, although Hatch et all 120051) have 
detected molecular line emission in several star-forming regions at 
distances of > 25 kpc from the center of the Perseus cluster. Our 
previous low resolution 4m A'— band spectroscopy has shown that 
there is a clear link between the intensities of emission lines from 
ionized and warm (~ 2000 K) molecular gas at the centres of these 
clusters. We now extend this work by exploring the kinematic and 
morphological link between these two phases to larger radii (and 
lower surface brightness limits). This is essential to separate the role 
of the AGN in these systems from that of the host galaxy and cool- 
ing flow in the excitation and dynamics of this gas. 

With the goal of connecting the cooler IGM states to that of 
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Abell 1795 


20 


21 03 2000 


3C445 


0.056 
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14.0 




9.1 




Abell 2204 


30 


29 07 2000 


Abell 3667 


0.056 


232 


<0.7 


<0.7 


<0.6 


<0.6 


Abell 2597 


40 


28 07 2000 


Sersic 159-03 


0.056 


233 


12.7 


12.7 


8.4 


8.3 


Sersic 159-03 


10 


13 07 2001 


Abell 2029 


0.076 


316 


< 2.9 


< 1.4 


< 4.8 


< 2.3 

* 








Abell 2204 


0.152 


632 


48.5 


* 


232.7 


Table 2. Observing Log 


of Clusters observed with Chandra 


Abell 2597 


0.082 


343 


41.5 


28.7 


80.7 


55.7 








Abell 1795 


0.063 


264 


*51.7 


N.A. 


*58.1 


NA. 









Table 1. Clusters observed with TTF, with observed total fluxes and lumi- 
nosities, or 3 <r upper limits. Fluxes are in units of 10~ 15 erg s _1 cm -2 . 
Luminosities in units of 10 40 erg s _1 . All values are uncorrected for ex- 
tinction. * No independent [Nil] measurement was taken for Abell 1795. 
The flux quoted for Ha includes the [Nil] emission. The [Nil] images of 
Abell 2204 were corrupted and could not be reduced. 



the hot X-ray gas on scales of > 10 kpc, we combine very deep 
narrow-band imaging of the BCGs of cooling flow clusters in the 
light of Ha A 6563 and [Nil] A 6584 lines with deep 8m long slit 
K-band spectroscopy of a subset, obtained in order to examine the 
distribution and thermal state of the molecular gas through the 1-0S 
and 2-1S emission lines of H2. The Paa emission in the K— band 
spectra allow us to link both data sets through comparison of the 
Ha and Paa lines. 



2 OBSERVATIONS 

2.1 Optical narrow-band imaging 

The Ha+[NII] observations were performed at the Anglo- Aus- 
tralian Telescope (AAT) using the Taurus Tunable Filter (TTF) 
to obtain images with a bandwidth of ~ 12A in all but 
one case. This Fabry-Perot e talon instrument, described in 
Bland-Hawt horn & Heath Jonej ll998h . has several advantages 
over conventional narrow band imaging: very narrow bandpass, 
very short wavelength step between "on" and "off bands and shuf- 
fling between on- and off- images on the CCD chip, which allows 
rapid switching between the two bands without increase in read- 
out noise. Together, these techniques yield images with very clean 
continuum subtraction, and independent maps of Ha and [Nil]. 

The galaxies observed are given in Table 1. These include 5 
strong cooling flow galaxies (as defined by the ROSAT observa- 
tions, before Chandra/XMM reobservations) and the non-cooling 
flow radio galaxy 3C445 (a known line-emitting source at a similar 
redshift to the clusters) and the non-cooling flow cluster A3667. 
Here and further below, apparent quantities are converted to in- 
trinsic quantities assuming H =72 km/s/Mpc, fl m = 0.3 and 
Q A = 0.7. 

All clusters were observed with a bandwidth of 12A (FWHP) 
except Abell 1795, which was observed with a bandwidth of 60A. 
Except for this cluster, each was observed at the redshifted wave- 
lengths of Ha and [Nil] (rest wavelength 6584A), although the 
[Nil] image for Abell 2204 was corrupted and could not be reduced. 
Abell 1795 was observed with a single broader bandpass covering 
both lines. For our data, the etalon step between on- and off-line 
images was equivalent to approximately 40A blueward of Ha and 
40A redward of [Nil] (about twice this value for Abell 1795 due 
to its wider bandpass). The off-line image was sufficiently far re- 
moved in wavelength to avoid picking up any extended line emis- 
sion present in any high velocity wing of the emission-line, but 



sufficiently close in wavelength that colour differences of the con- 
tinuum between on- and off- bands will have a negligible effect on 
its removal from the on-line image. 

As expected from previously-published less-sensitive obser- 
vations, optical emission lines were detected from five of the 
sources. No significant emission was detected from Abell 2029 
or the control cluster Abell 3667. Grey scale and contour plots 
of the detected clusters are shown the left and center panels of 
Figures 1 1151 respectively and for four clusters the contour plots 
are overlayed on Chandra X-ray images in the righthand panel. 
These images were extracted from the Chandra science archive 
site http://asc.harvard.edu/targets They represent data from the 
ACIS-S3 chip, summed over the 0.3-7.0 keV energy range, and re- 
duced with the CI AO 3 . software package. The photon counts 
were smoothed with a 4.4 arcsec full-width at half maximum gaus- 
sian kernel. The log of the Chandra observations is given in Table 2. 

In Table 1 we also list the integrated fluxes and luminosities 
contained within a perimeter drawn at a level equal to the rms noise 
in the image. In Figures [6] and [7] we plot the azmuthally averaged 
surface brightness of the emission lines in each of the four well 
resolved sources as a function of radius from the galactic nucleus. 



2.2 A' band spectroscopy 

In order to probe in detail the connection between the atomic and 
molecular phases of the line-emitting gas, we obtained medium- 
resolution spectroscopy of a subset of the above cluster sample. 
These observations were obtained with ISAAC on the Very Large 
Telescope (VLT). The sample was that subset of the optical imaging 
sample described above, showing emission at large radii and easily 
visible from the Southern hemisphere. FiguresQ[2]and|4]show the 
positions of the slits superimposed on the contour maps obtained in 
the previous section. 

Because of the extended nature of the Ha emission we used 
a relatively wide slit (1") for the IR spectroscopy, assuming (cor- 
rectly) that the IR lines would be similarly extended. The slit was 
generally oriented through the galaxy nucleus and at a position an- 
gle at which we had previously observed extended Ha emission. 
We used the ISAAC medium resolution (MR) grating, with a reso- 
lution of ~ 3000 after experience showed that the results with the 
low resolution (LR) grating (R ~ 500) were inferior. The main ad- 
vantage of the MR system is that the K-band contains several strong 
telluric absorption systems near the 2-micron end of the band and 
with the LR system the individual lines in these systems blend to- 
gether. This makes it difficult to detect weak extraterrestrial emis- 
sion in this region. With the MR system the individual lines are 
separated, and emission can easily be detected between the absorp- 
tion lines. 

A disadvantage of this approach is that we cannot observe 
the entire K-band in one setting of the spectroscope. The observ- 
able band was about 0.12/i at any one setting of the central wave- 
length. We have chosen to take one spectrum in the vicinity of Paa, 
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Figure 1. Images of the line emission around A2597. North is up and East at left. Coordinates give the offsets from the nucleus in arcsec on the sky. At 
Abell 2597, 1 arcsec = 1.66 kpc. The small feature in the bottom right is the residual image of a bright star. Left: Greyscale and Middle: contour plot of 
Ho+[NII]. Lowest solid contour is at 3. Ox 10~ 18 ergs s _1 cm _2 arcsec -2 . In this and following figures, the subsequent contours double in surface brightness, 
while dashed contours represent negative (noise) values. Right: Contour plot of line emission overlayed on greyscale of Chandra X-ray emission from cluster. 
For clarity in the overlay every other contour has been omitted. 






Figure 2. Images of the line emission around Sersic 159-03. Left: Greyscale and Middle: contour plot of Ha+[NII]. Lowest contour is at 3.0x10 18 ergs 
s _1 cm _2 arcsec~ 2 . Right: Contour plot overlayed on Chandra image. 1 arcsec = 1.13 kpc 





Figure 5. Contour plot of 3C445 H« emission. Lowest contour is at 5.25 X 10 18 ergs s 'cm 2 arcsec 2 . 1 arcsec = 1.13 kpc 



H2:1-0S(3) (and occasionally other lines), and one in the vicinity 
of H2:1-0S(1). Because of variations of the exact bandpass with 
redshift, and small errors in redshift, some of the less important 
lines at the edges of the spectra (e.g. Br7 were not observed for all 
clusters). 

Standard stars, usually B-stars of magnitude K~ 7, were ob- 
served to correct for telluric absorption and to provide calibratable 
fluxes. These standard spectra were used to remove the continuum 
emission by a linear fit to the galactic emission in regions where 
no galactic emission or absorption lines are expected. Additional 
smooth power-law components are also removed, which mostly re- 
flect the difference in blackbody slope between the galactic stars 
(mostly K- and M- giants) and the standards. Any absorption lines 
in the standard stars then remain as emission lines in the final spec- 
tra. This is sometimes a problem near rest wavelength Br7. The 
aforementioned spectral mismatches appear mostly as very weak 
absorption or emission features, before correcting for the total tel- 
luric absorption. 

These "raw" two-dimensional spectra are shown in Figures|8| 
through [TOl In these 2d spectra we see the gaseous emission lines, 
along with noisy vertical stripes at wavelengths of high telluric ab- 
sorption. By Kirchoff 's principle, the atmospheric thermal emis- 
sion at these wavelengths is maximum. The photon noise from the 
atmospheric emission is the chief source of noise in the final spec- 
tra, and thus maximum at the position of the absorption lines. At the 



y-position of the maximum of the galactic continuum we also see 
small features due to minor stellar absorption/emission mismatches 
between target galaxy and standards. 



3 DISCUSSION 

3.1 Extent, morphology and kinematics of the "warm" gas 
phases 

Figures Q through [5] show that the warm HII phase (T ~ 10, 000 
K) can be traced to radii beyond 20 kpc in all four cooling flows. In 
Abell 2204 we see emission to 30 kpc and in Abell 2597 to beyond 
40 kpc in smoothed images. This last figure can be compared to HII 
structu re seen in the same cluster to about 20 kpc in He ckman et alJ 
119981) and the bright central features associated with the c entral 
radio source at ~ 5 kpc shown in Koekemoer et a D 1 19991) . The 
increased maximum radius of detection for this source corresponds 
to an incr ease in surface bright ness sensitivity by a factor of 20 
relative to lHeckman et alJ ll998t) . 

The four clusters have diverse morphologies: extended and 
roughly circularly symmetric (Abell 2597), compact (Abell 2204), 
and arc-like (Abell 1795 and Sersic 159-03). Ignoring these dif- 
ferences, we have plotted the radial surface brightness distribution, 
averaged over annuli centered on the galaxy nucleus in figures |5| 



HII and H 2 in cooling flows 5 



Abell 2597 Ha + [MI] 



Sersic 159-03 H«+[NH] 



Power Low (-3) 
Power Low (-2) 



■II. 



Power Law (-2) 
Power Law (-3) 



'fi. 



T . T 



Figure 6. Left: Averaged surface brightness in combined Ha and [Nil] emission in circular annuli around the nucleus of Abell 2597. Right: Sersic 159-03. 
The error bars include noise, intrinsic non-circularity, and estimates of the uncertainty in the background subtraction at large radii. For radii less than ~ 15 
kpc the eiTor bars are smaller than the plot symbols. Also shown are power law curves for r~ 2 and r~ 3 laws and for Sersic 159-03 an exponential curve 
£ oc exp(— r/a). 



Abell 2204 Ha Abell 1795 Ha + [NII] 




Figure 7. Averaged surface brightness in circular annuli around the nucleus of Left: Ha in Abell 2204; Right: Ha combined with [Nil] in Abell 1795. Also 
shown are power laws, for Abell 2204 a gaussian relation: £ oc exp(— r/a) 2 , and for Abell 1795 an exponential relation. 



and [7] The averaging over a large number of pixels at larger radii 
improves the signal/noise ratio at these radii. 

Abell 2597 shows significant emission, distributed essentially 
as an r~ 3 power law to beyond 40 kpc, where it disappears 
smoothly into the noise. The "arc-like" structure of Abell 1795, 
when azmuthally averaged, shows similar behaviour, with perhaps 
a slightly flatter distribution with radius (oc r~ ). Abell 2204 shows 
emission to a similarly large radius (~ 35 kpc) but drops more 
rapidly than a power law, more like a gaussian, with a radius of 
~ 12 kpc. In other words the apparent size of the emission re- 
gion seems to be a real physical characteristic of the source, rather 
than a limit imposed by our surface brightness sensitivity. In Ser- 
sic 159-03 the circularly averaged emission can be traced to radii 
just beyond 20 kpc and is fit approximately by an exponential with 
scale sizes of ~ 4 kpc. Thus the increased sensitivity afforded by 
the TTF system yields two contrasting results: in Abell 2597 and 
Abell 1795 we can trace the smooth decrease in surface brightness 
to double the previous limiting radius but find no physical edge; 
in the other two clusters we can establish the existence of a radius 
characteristic of the source physics. 

Figures ISl 1 01 show strikingly that molecular emission can be 
traced to ~ 25 kpc and that the Paa and H2 line profiles are virtu- 
ally indistinguishable both spatially and kinematic ally. A compari- 
son of the Paa lines and the 1-0S(1) and 1-0S(3) lines shows that, 
taking differing signal/noise ratios and interfering sky lines into ac- 
count, every feature visible in one line is visible in the others at 
the same velocity and position and essentially the same intensity. 
This implies directly that that the HII and warm H2 gas phases 
and their respective excitation sources must be well mixed on scale 
sizes finer than 1 kpc, and the relative mix is essentially constant 
over factors of 10 in radius. The correspondence suggests further 



that the H2 phase is present wherever the diffuse HII can be ob- 
served, i.e. to beyond 40 kpc in Abell 2597. 

Consistent with our earlier work, the ratio of H2 to HII line 
strengths is far higher than that found in typical galactic star form- 
ing regions and starburst galaxies (by HII lines we mean hydrogen 
recombination lines from the Balmer and Paschen series, arising in 
HII regions). In fact the combined luminosity in all the Near In- 
fraRed (NIR) H2 lines is comparable to that in all the optical HII 
lines, although the luminosity in Lytv is probably an order of mag- 
nitude above this. 

The velocity structure of the molecular gas in all three ob- 
served clusters shows a division into two regions: a central area of 
radius ~ 5 kpc with high velocity dispersion and a larger region of 
5-25 kpc with much lower dispersion. The one-dimensional spec- 
tra presented in previous works including our own are dominated 
by emission from this central region. 

In this region, where the gas is probably stirred by the AGN, 
the velocity spreads - here quoted as the line-of-sight rms velocity 
dispersion a - range from ~ 120 km s" 1 for Sersic 159-03 through 
200 km s" 1 for Abell 2597 to 300 km s" 1 for Abell 2204. The ar- 
eas at larger radii show lower dispersion (< 60 km s _1 ) and appar- 
ent large scale rotation or shear motions of magnitude ~ 100 km 
s" 1 . These dispersion estimates are in fact upper limits set by the 
velocity resolution of our wide slit spectroscopy and the true val- 
ues could be lower. For Abell 2204 and Abell 2597. lEdgeN200ll) 
finds line widths in the CO molecular emission of 255 and 278 km 
s _1 respectively with a beam size of 22 arcsec. These values seem 
lower than we would expect from our K-band data, given that our 
central FWHP "widths" are 700 and 470 km s _1 for these two clus- 
ters, and that the K-band line emission is strongly peaked toward 
the centres. This would indicate that the cold molecular gas respon- 
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Figure 8. Raw two-dimensional spectra of Abell 2597 in the regions of Pao and 1-0S(1) respectively. The abscissa gives the laboratory wavelength in microns. 



sible for the CO emission does not partake of the violent motions 
near the galaxy center and is situated at radii beyond ~ 10 kpc. 

At the larger radii the gas velocity dispersions and rotational 
velocities are much smaller than the typical stellar dispersions ve- 
locities in BCGs: 400-600 km s" 1 . Thus the the molecular gas 
clouds at large radii are not kinematically supported nor, given their 
high density, can they be supported by pressure gradients in the sur- 
rounding hot gas. The low velocities that we see are dynamically 
consistent only if this gas has been formed from a pressure sup- 
ported phase in the last 10 7 years (the time necessary for gravity to 
accelerate them to abo ve 100 km s -1 ) or an as yet unknown mech- 
anism supports them. lDaines et al]ll994l) have suggested magnetic 
support of the clouds. This seems unlikely to us however, because 



of the Parker instability IParkeill 19681) generally present in gas sup- 
ported against gravity by magnetic fields. 

From our Ha fluxes, with corrections for dust absorption (see 
below) we can calculate the total mass of warm ionized gas in the 
clusters centres. Because of the relatively fast drop of brightness 
with radius we do not, however, expect these estimates to be much 
larger than older data obtained at higher surface brightness limits. 

F or Abell 2597, assuming n e ~ 200 cm -3 tHeck man et all 
1998) and A Ha - 1 we find an HII mass of ~ 3 x 10 7 
M©, The mass in warm H2 can be estimated from Mh 2 = 
2m,HL\/(eiA\Fi) where L\ is the luminosity in the 1-0S(1) line, 
ei is the energy of the photon from this transition, A\ is the Ein- 
stein coefficient for this transition, and Fi is the fraction of all H2 
molecules to be found in the upper state for this line. For LTE con- 
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Figure 9. Raw two-dimensional spectra of Sersic 159-03 in the regions of Pact and 1-0S(1) respectively. 



ditions at T ~ 2000 K (c.f. Section l3!24l . Fi is ~ .015. With the 
assumption that the ratios of 1-0S(1) to Paa and Paa to Ha vary 
globally the same way as they do along our slit measurements, we 
find Li ~ 10 42 erg s' 1 for both Abell 2597 and Abell 2204, with 
a corresponding Mh 2 of ~ 4 1 s Mp. This is much smaller than 
the masses of cold H 2 found bv lEdeel EooH) to be ~ 10 9 " 10 M . 
The H2 masses calculated in that a rticle are based on "standard" 
ratios of H2/CO and llsraell Il997l) shows that these may in fact 
be serious underestimates in cases where the molecular gas is ex- 
posed to a bright UV/X-ray radiation field, which dissociates the 
CO molecule. The kinetic energy in the warm molecular gas in the 
"violent" central region is of order 10 54 ergs and if it is dissipated 
in the central dynamic time of ~ 2 x 10 r years the corresponding 
kinetic luminosity is of order 10 39 erg s -1 , which is small relative 



to other luminosities associated with the AGN and cooling flow. If, 
contrary to our inference above, the cold gas also partakes of this 
rapid turbulent motion, its kinetic luminosity would of course be 
much larger and would become important to the energetics of the 
system as a whole. 

The problem of the support of the molecular gas takes on its 
most serious form for the cold "CO" gas. The gravitational poten- 
tial energy of this gas,<E> co zd can be estimated as M co id<ri, where 
a, is the 3-dimensional velocity dispersion of the stars ~ 600 km 
s 1 ; 4>coid ~ 3 x 10 58 erg. This is the same order of magnitude 
as the total th ermal en ergy of the X-ray gas within 15 kpc of the 
galaxy center jKaastra et alJl2004l) and similarly the entire mag- 
netic e nergy in this region, assuming B ~ 10/iG (e.g. Tayl or et alJ 
(2002])). 
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Figure 10. Raw two-dimensional spectra of Abell 2204 in the regions of Paa and 1-0S(1) respectively. 



3.2 Line ratios and thermal structure of the warm gas phases 

Figures l8l to 1 1 31 illustrate the similarity of the surface brightness 
profiles for the NIR HII and H2 lines. To allow a more direct com- 
parison of the line ratios of the strongest lines (1-0S(1), 1-0S(3) 
and Paa, we have integrated them in velocity and plotted them as 
a function of slit position in the third panels of Figures H31to[Tg| In 
each case the spectra have been plotted on a logarithmic intensity 
scale, and scaled to the flux of 1-0S(1) at the center. In the fourth 
(d) panel of these figures we have plotted in a similar fashion the 
intensities of Ha, [Nil] and Paa , where the optical lines have been 
evaluated from the TTF images along the position of the NIR slit. 
In estimating the absolute IR/optical ratios, some uncertainties in 
this procedure are inevitable: the different observing and calibra- 



tion techniques may introduce errors of ~ 30%, and errors in the 
exactly alignment of the slit on the images may be of order 1 arcsec. 
Experiments with moving the positions or position angles of the slit 
on the image did not produce significant changes in the shapes of 
these plots. Additionally, the narrow bandpass (12A) of the TTF 
filters will miss some optical flux at the central high velocity width 
position. Assuming a gaussian profile with the dispersions plotted 
in the Figures, we find about 40% of the Ha flux at the center of 
Abell 2204 would be lost (cr=300 km s" 1 ), 20% at the center of 
Abell 2597 (cr=200km s" 1 ), and a negligible loss in Sersic 159-03. 
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Figure 11. Contour maps of calibrated stronger detected lines in Abell 2597. Contours are [-1,1,2,4,8,16...] times the lowest level (in erg s _1 cm — 2 arcsec - 2 ), 
which is indicated in the figure title. The vertical bar in the 1-0S(1) map marks the regions for which LTE line ratio plots are given in FigurefT8l 



3.2.1 Dust extinction 

Under the assumption of Case-B ionization equilibrium in the HII 
regions, the Pq/Hq ratio is ~ 0.12 l Osterbrock 1974). All three of 
the objects with ISAAC spectra show larger ratios at the peak of 
their surface brightness profiles: 0.76 for Abell 2597, 0.55 for Ser- 
sic 159-03 and 0.24 for Abell 2204. For the first two clusters this 
implies extinction of Ana ~ 1.6 — 1.8 mag with calibration un- 
certainties of ~ 0.5 mag, and reddening of E(B — V) ~ 0.6. For 
Abell 2204, after correction for Ha lost outside the filter bandpass, 
the estimated extinction is small, < 1 mag. In determining these 
numbers, we have assumed that the dust is in a screen in front of 
the gas, whereas in reality the dust could well be mixed-in. Never- 
theless,the numbers give an indication of the amount of dust in or 
around the emission-line regions. 

In the case of Abell 2597 this ratio decreases smoothly by a 
factor of 2 — 3 as r increases to ±15 kpc from the centre of the 
nebulosity, corresponding to a drop in extinction to ~ 1 mag. We 
see a similar trend for Abell 2204 probably due to a reduction in 



loss outside the Ha filter i.e. showing little evidence for extinction 
at any radius. 

The situation is more complicated for Sersic 159-03. Along 
the "tail" of emission north of the peak surface brightness, the line 
ratio barely varies significantly, possibly the Pa/Ha ratio increases 
slightly relative to the peak, implying an increase of 0.2-0.3 in E(B- 
V) at most. To the south, the surface brightness drops sharply. Al- 
though Figure fT5h shows a sharper drop in Paa than Ha indicating 
a decrease in reddening, differences in the seeing and pixel size be- 
tween the two data sets may cause this. These systematic effects 
limit our ability to draw firm conclusions about the effect of red- 
dening in this region. 

What is clear from the above is that many line emitting clouds 
up to 20 kpc of the centres of these clusters contain dust and the red- 
dening effect of the dust does not vary dramatically across the neb- 
ulosity, apart from at their very centres. iKoekemoer et alj 1 19991) 
have already shown that there is a probable dust lane crossing the 
nucleus of Abell 2597, this is the obvious explanation for the excess 
reddening we see at the centre for this source. 
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Figure 12. Contour maps of calibrated stronger detected lines in Sersic 159-03 
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Figure 13. Contour maps of calibrated stronger detected lines in Abell 2204 



3.2.2 [NIl]/Ha line ratios 

The ratio of the forbidden [Nil] line to the permitted Ha line 
can be taken as a crude measure of ionization parameter H = 
n-ohotonln ii; larger values of rNI Il/Hq correspond to larger val- 
ues of H. Hec kman et alj 1 19981) have shown that from clus- 
ter to cluster, various ionization measures correlate with each 
other, and that ionization correla t es inversely with Ha lumi- 
nosity. ICrawford & Fa bian ( 1993), ICrawford et alj J 19951) and 
ICrawford et alJ 1 19991) have shown that these correlations are 
quite general in cooling flows, but that the line ratios vary con- 
tinuously rather than forming two distinct classes as claimed by 
iHeckman et alj (l998). For the two clusters where we have good 
TTF maps of these lines, the line ratios are displayed in Figure [TTl 
We find that between these two clusters, the relation between 
[NII]/Ha varies in the sense described in the quoted papers, but that 
locally within each cluster, this is not true, and that the two clus- 
ters appear quite different. In Abell 2597 there is a some tendency 



for this ratio to increase where the surface brightness of Ha de- 
creases, but the effect is only strong at the very edges of the nebulae. 
Log([NII]/Ha) is in fact very nearly constan t at a low value ~ — 0.2 
typical of the high luminosity "Type II" of IHeckman et alj \ 19981) 
along a NE-SW axis, more or less parallel with the radio jet axis, 
rising to the "Type I" value of ~ +0.3 on the NW and SE perime- 
ters. The source average is dominated by the central, lower value. 
Sersic 159-03 shows the opposite behaviour, the brighter central 
regions show a high value of log([NII]/Ha), while the perimeter 
shows lower values; the high value dominates the average. 

We note that the ionization level is not simply related to dis- 
tance from the galactic nucleus. In Abell 2597 the line ratio varies 
little when the radius varies by a factor of more than 10. This fact 
emphasizes, as other authors have noted, that the source of excita- 
tion is almost certainly not the AGN but rather of an origin local 
to the ionized gas. The difference between the two clusters empha- 



HII and H 2 in cooling flows 1 1 



A 
> 
V 



200 
150 
100 
50 


-50 



Abell 2597 1 -0S(3) 



200 

150 

100 

50 
0. 



< 

E 



cn . „-i7 



i ; i 

+ 1-0S(1) 

O Paa * 2.74 

* 1-0S(3) * 1.02 




# <> 



K ' 0" 



D 



+ Ha 

^ Paa* 1 .32 



iff; 




< 

> 


> ° 











-30 



-20 



■10 




kpc 



20 



JJ 



30 



Figure 14. Abell 2597 spectral data as a function of distance along the slit: top and second panels give mean velocity and velocity dispersion taken from the 
1-0S(3) line. The third panel gives the surface brightness of the IR lines. The ordinate gives the surface brightness of the 1-0S(1) line. The other lines (Paa 
and 1-0S(3)) have been scaled by the amounts indicated in the figure to match l-OS(l) at the center. The points for all three lines are measured at the same 
positions, but the Paa points have been shifted 0.5 kpc to the right, and 1-0S(3) 0.5 kpc to the left for clarity. The bottom panel gives the surface brightness of 
ionized lines. The ordinate gives the surface brightness of Ho while Paa has been scaled to match Ha at the center. 
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Figure 15. Sersic 159-03 spectral data: top and second panels give mean velocity and velocity dispersion taken from the l-OS(l) line.The third panel gives 
the IR lines; the ordinate gives the surface brightness of the l-OS(l) line and the other lines have been scaled by the amounts indicated to match l-OS(l) at 
the center. The points for the various lines have been shifted slightly in the x-direction as described in the previous plot. The bottom panel gives Ha and Pact 
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Figure 16. Abell 2204 spectral data: top and second panels give mean velocity and velocity dispersion taken from the Pa alpha line. The third panel gives 
1-0S(1) and Paa scaled to match 1-0S(1) at the center. The bottom pannel gives Ha and Paa scaled to match Ha at the center. 



sizes that the interpretation of these line ratios is not simple, as they 
depend on metallicity, ionization parameter, and density. 



3.2.3 Comparison between H 2 lines and Paa 

Using our ISAAC spectroscopy, we are able to compare the spa- 
tially resolved properties of the Pa and H2 1-0 lines. Previous work 
has concentrated on the integrated properties of the highest sur- 
face brightness IR lines. These spectra are invariably dominated 



14 Jaffe, W., Bremer, M.N., Baker, K. 





Figure 17. (left) Greyscale representation of log([NII]/Ha in the central region of Abell 2597 with superimposed the contours from FigureQ The orientation 
is the same as FigureQ while the field of view covers 60" in both coordinates. Where the signal/noise was low the greyscale has been blanked to a uniform 
grey, (right) log([NII]/Hct for Sersic 159-03 with the contours from Figure|2] Again a 60" field is plotted. 



by emission from the very centres of the galaxies (which harbour 
strong radio sources) and likely to be influenced by the effect of 
AGN emission. 

We can compare three main spatially-varying parameters of 
the lines: Their surface brightness profiles, their velocity disper- 
sions and their peculiar velocities relative to the centre of the 
galaxy. Figures 6-14 show that the atomic and molecular hydrogen 
lines behave very similarly as a function of position. 

The velocity structure in all the lines looks similar. As noted in 
Section l3~Tl fhe central high surface brightness regions have veloc- 
ity dispersions of several hundred km s" 1 , while at larger radii the 
dispersion in all lines drops to approximately 60 km s -1 , the reso- 
lution of the spectrograph, and probably lower. The molecular and 
atomic lines show similar variations in peculiar velocity as a func- 
tion of position, showing that these gas phases are closely linked 
kinematically, even though their dynamical properties, e.g. the ra- 
tio of gravitational to drag forces, are quite different because of the 
large difference in den sity: n e < 10 3 cm -3 in the ionized phase 
Heck man et alj ^,998) versus n e ~> 10 5 cm~ 3 in the molecu- 
lar phase (c.f. Section [3.2.41 . The surface brightness profiles of 
the atomic and molecular phases also follow broadly similar pat- 
terns. In Figures fT4b to ll 6b the brightness profiles integrated across 
the slit are normalised to the brightest point in the 1-0S(1) profile. 
For Abell 2597 the ratio Paa/1-0S(1) appears a factor two lower at 
the peak of the profile than in the rest, with molecular and atomic 
emission otherwise tracing each other. For Abell 2204 the ratio de- 
creases by no more than a factor of two from the centre to the edge 
of the detected nebulosity. For Sersic 159-03 the ratio appears con- 
stant across the northern tail, with a decrease of a factor of 2 in the 
relative strength of Pa at the peak of the profile, similar to that seen 
in Abell 2597. The current data do not allow us to trace the atomic 
emission to the south, but the non-detection is consistent with the 
ratio found for the north. 

We have earlier noted Ijaffe et all 1200 ll) that this ion- 
ized/molecular line ratio is relatively constant from cluster to clus- 
ter, and is typically at least a factor of ten lower than that usually 



seen in photo-dissociation regions (PDRs) around star-forming re- 
gions. The current observations show that this ratio does in fact vary 
somewhat, by factors of ~ 2 from cluster to cluster and within clus- 
ters. We will discuss in Section l3~3l some of the possible reasons for 
this variation, but our principle remark is that the excitation sources 
of the HII and warm H2 phases must be very closely linked in order 
that the line ratios and velocities remain so close to each other de- 
spite changes in radius factors > 10, surface brightness by factors 
> 100 and local velocity dispersion by factors > 4. 

If we accept the hypothesis that the ratios of HII and H2 lines 
is essentially constant over the area that we see HII line emis- 
sion, we can calculate total galaxy luminosities in the H2 lines. As 
given in SectionO for the 1-0S(1) line in Abell 2597 this yields 
Li ~ 2 x 10 42 erg s~\ for Sersic 159-03 Li ~ 2 x 10 41 erg 
s _1 and for Abell 2204 L\ ~ 4 x 10 41 erg s _1 with uncertainties 
of order a factor of 2. The luminosities in the H2 1-0S(3) are sim- 
ilar, and detailed mod els of excitation of H2 gas in many c ircum- 
stanc es (e.g. [Black & van Dishoecllll987tlsternberg & Dalgarno. 
1989) generally find that the luminosity in the strongest single ro- 
tational/vibrational line is generally ~ 1 — 2% of the total molec- 
ular NIR emission due to H2. This implies a typical total H2 NIR 
luminosity of 10 43 — 10 44 erg s _1 . This is up to two orders of mag- 
nitude greater than the Ha luminosity and may be greater than the 
total optical line luminosity (assuming that Lya luminosity is ~ 10 
times that of Ha). 



iKaastra et aT]j2004 quote total X-ray fluxes (0.1-2.4 keV) for 
Sersic 159-03 and Abell 1795 that are approximately 2000 times 
larger than the Ha fluxes quoted in Table 1. With the conversion 
factors used in the previous paragraph the total NIR luminosity for 
these galaxies is ~ 5% of the total X-ray luminosity at all radii, and 
a larger fraction of that arising in the central ~ 20 kpc. For example 
in the Chandra data of Abell 1795 jEttori et alj|2002l) only about 
4% of the total X-ray emission in the 0.5-7 keV range arises in the 
central 17 kpc. This suggests that the NIR and X-ray luminosities 
in this central region are similar. 
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Figure 18. LTE diagrams for Abell 2597. Different regions along the spectrometer slit, marked in Figure HTI are given in different panels above. The abscissa 
marks the excitation energy of the upper state of each transition (expressed in temperature units), while the ordinate gives the flux of the corresponding line, 
divided by its Einstein coefficient, A and the statistical weight of the upper state g. Finally each flux is normalized to the flux of the 1-0S(1) line. In LTE, 
F oc gAexp(—Eu/kT), so the various transitions should lie along a straight line. 
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Figure 19. LTE diagrams for Sersic 159-03. The 1-0S(4) transition lies very close to a sky line at this redshift, and is not given in the plot. 



3.2.4 H2 line ratios and LTE calculations 

In FigureQUwe show the relative strengths of the 1-0S(1),S(3),S(4) 
and 2-lS(3) lines averaged within five spatial regions marked on 
Figure fTTI and displayed in the form of log(F/gA) versus E upper 
where F is the line flux, g the statistical weight of the upper state, 
A the Einstein coefficient of the upper state, and E upper the energy 
of the upper state of the transition expressed in temperature units. 
In local thermodynamic equilibrium (LTE), the resulting points 
should lie on a straight line with slope 1/Tlte- For Sersic 159- 
03 the same plots are given in Figure [T9l but the 1-0S(4) line is not 
included because its wavelength falls near that of a bright sky line. 
The values and uncertainties of Tlte printed in each plot are de- 
rived solely from the 1-0S(1) and S(3) lines, which have by far the 
best signal/noise ratios. 

Looking first only at the temperature estimates from these two 
lines we see in Abell 2597 a gradient of decreasing temperature 
from the the southern region 1 to region 4 at a height of ~ 10 
kpc. While it is tempting to conclude that the temperature decreases 
with distance from the nucleus the estimated average value for the 
southern region 1 is actually higher than that for region 2. The tem- 
perature uncertainty in region 1 is, however, quite large due to the 
low signal level there. Comparison of the error bounds indicates 
that while T\ could be the same as T2, it cannot be not significantly 
lower, e.g. not as low as T3. The weighted best LTE temperature 
for all five regions together is ~ 2600 K, somewhat hi gher than the 
central value of 2160 K given bv lWilman etai]fc002l) for this clus- 



ter, but within the range typical for cooling flows. In Sersic 159-03 
there is an indication of increasing temperature with distance from 
the nucleus, but the significance is low. The global average temper- 
ature along the slit is ~ 2400 K. 

In none of the regions is there significant evidence for line ra- 
tio deviations from the LTE curve. There are several regions where 
the signal at the 2-lS(3) line, which has an excitation energy of 
13890 K, is significant, namely regions 2 and 3 of Abell 2597 and 
region 1 and 2 of Sersic 159-03. In these cases the positions of this 
line near the LTE line imply densities of n ~ 10 7 cm -3 in order 
for there to be enough energetic collisions to keep the upper level 
populated. 

3.3 Source of molecular excitation 

The problem of the source of excitation of the ionized gas, and 
to a lesser extent that of the molecular gas, has previously been 
discussed by many authors, including ourselves. The new contribu- 
tions of the data reported here consist firstly of the demonstration 
of the very close spatial and kinematic correlation of the HII and 
H2 lines and secondly of the extension of these correlations to large 
radii with respect to the central galaxy. 

The strong correlation between the emission characteristics of 
the two gas phases argues that they are excited by a single physical 
mechanism and not, for example, that stellar photoionization ex- 
cites the HII region, while heat conduction from the X-ray medium 
excites the H2 emission. Shock excitation provides such a natu- 
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ral connection because the various gas phases represent a time se- 
quence in the shock heating and cooling hist ory, but opt ical spectra 
in cooling flows at t he rOIIIl A 43 63 line IVoit & Donahuel ll99^ 
and more stringently Baker 1 2004)) are not consistent with shock 
models, and we will not discuss these further. 

A more likely explanation is excitation of both phases by 
relatively hard radiation. In starforming photodissociation regions 
(PDRs) excitation by young stars causes HII line emission due 
to absorption of Lyman continuum photons and subsequent re- 
combination. In these regions, H2 lines arise from the absorption 
of Werner and Lyman lines by H2 at longer wavelengths near 
1000 A and radiative or collisional de-excitation of the mol ecules 
(Blac k & van DishoecklfT987llSternberg & Dalgarnclll989l) . Pho- 
tons beyond the Lyman limit have a low cross section for absorp- 
tion in the HII regions, but can be absorbed in dense H2 regions 
with the release of energetic photoelectrons which lo se their en- 
ergy primarily by spinning up H2 molecules (Tineetal., 1997; 
Malo nev et all 1 1996). This efficient conversion of hard photons 
to molecular vibration-rotation lines should lead to higher molecu- 
lar/ionic line ratios than those from softer sources. In fact, in H2 re- 
gions surrounding planetary nebulae, the molecular/ionic line ratios 
are similar to those in cooling flows (Jaffe and Mellema, in prepa- 
ration). We have not yet carried out detailed modelling of these line 
ratios in cooling flows, because the large pressure difference be- 
tween the two phases implies that a full dynamic model, including 
outflows is probably necessary. 

The assumption of a hard ionizing spectrum is attractive be- 
cause it is consis t ent with detailed analyses of the optical spectra 
IVoit & Donahue! Il997t Erawford & Fabiail Il992t iBakeJ. |2004 
which indicate the presence of radiation from a source harder 
or hotter than O-s tars (e.g. white dwarfs or Wolf-Rayet stars). 
lO'Dea et alJ 120041) find evidence of copious FUV radiation in the 
regions of Abell 1795 and Abell 2597 that show molecular emis- 
sion in our maps, although the radiation detected in the HST/STIS 
observations is not "hard" compared to Ly man continuum. The 
source of these hard photons is still unknown. Crawford & Fabian 
1 1992) suggest, for example, radiation from a non-equilibrium 
"mixing-layer" of X-ray and molecular gas. Our observations re- 
enforce the need for hard radiation, but do not give direct evi- 
dence for its source. They do however lead to several interesting 
constraints on the spectrum and spatial distribution of the exciting 
photons. 

The cooling rate of dense H2 gas at LTE from rota- 
tional/vibrational radiation alone is qui te high and the gas is quite 
difficult to keep warm (e.g. ICrawford & Fabianl fl992). The cool- 
ing rate per molecule from this radiation is about 3 x 10~ 17 erg 
s _1 near 2000 K. The heating rate is Fecth where Fe is the ex- 
citing flux (erg s _1 cm -2 ) at energy E, and the cross sec- 
tion for absorbing photons of energy E, which is roughly ah ~ 
5 x 10~ 18 (S/20eV)" 3 cm 2 . If we assume 100% efficiency of con- 
verting hard photons to NIR emission, and a uniform distribution 
of the heating sources within a radius r ~ 20 kpc, we can estimate 
F E ~ L NIR /A-Kr 2 . For L NIR ~ 3 x 10 43 erg s"\ F E ~ 10 -3 
erg s _1 cm -2 . Thus even for photons with energies just above 1 
Rydberg (say 20 eV) the heating rate per molecule is ~ 5 x 10 -21 
erg s _1 or four orders of magnitude below the cooling rate. 

Thus if photons beyond the Lyman limit are the heating 
source, we must conclude both that these photons are relatively soft 
(EUV rather than soft X-ray) and that the surface filling factor of 
the photon sources is very low. In other words, the sources of pho- 
tons must be clumpy and lie very close to the molecular gas, so that 
the local value of Fe is much larger than its average value. Be- 



cause of the difficulty of detecting EUV photons, there is currently 
little direct evidence fo r such emission in cl usters and existing ob- 
servations are limited. TOegerle et alJ (2001) detected O VI emis- 
sion from A2597 in FUSE observations, indicating the presence of 
gas at ~ 10 5 K., whereas fLecavelier des Etangs et alJ 12004) failed 
to detect the same lines in Abell 1795. Although a soft excess of 
EUV and soft X-rays has now been detected in several clusters 
(e.g. QCaastra et al. 2003J), this is generally detected across the en- 
tire cluster and only in the case of Sersic 159-03 shows any sign 
of peaking close to the cluster centre. The large-scale of the emis- 
sion contrasts with the highly clumped nature of any EUV heating 
source indicated by our data. 

The CO observations of lEdgd fcOOll) support the idea of 
clumpy, low surface filling factor, molecular gas; for several clus- 
ters the absorbing column calculated from an assumed smooth dis- 
tribution exceeds that found in X-ray observations. For th e clusters 
Abell 1068 and Abell 1835, for example. Edge] 1200 ll) find val- 
ues of Iqo d-o) = f T a dV of approximately 3 and 2 K km s _1 
respectively. Ilsraell ll997l) shows that the ratio of the column den- 
sity Nh 2 (cm -2 ) to ico(l-o) i s typically ~ 3 x 10 21 in external 
galaxies. Thus the CO emission fluxes are equivalent to column 
densities Nh 2 of 9 and 6 xlO 21 cm -2 . The estimated X-ray ab- 
sorption column densities are , in contrast, about 2 x 10 20 for both 
clusters JWise et all2004l and lPeterson et all200ll) . The molecular 
gas implied by the CO observations would have been visible in the 
X-ray observations had it smoothly covered the regions of X-ray 
emission; the absence of this absorption implies a small covering 
factor for the molecular gas. 

The luminosity needed to heat the H2 gas an d that missing 
from the relatively cool gas in cooling flow models IPeterson et alJ 
120031) are both of order ~ 10% of the total central X-ray flux. This 
prompts us to speculate that the H2 medium in fact absorbs and 
hides most of the softer radiation from the cooling flow. While the 
X-ray spectra are not consistent with a foreground screen of absorb- 
ing gas covering all the emitting gas, they might still be consistent 
with clumpy selective absorption. Here we postulate that the cooler 
region s of the X-ray gas, with T < T am i, ien t/3 IPeterson et alJ 
2003) consist of small clouds each surrounded by even colder, 
dense, H2 gas. Because many other processes in the central regions 
of rich clusters have similar radiative or kinematic luminosities, 
however, the issue of the H2 excitation source and/or the heating 
of the cooler X-ray gas is still wide open. 



4 SUMMARY AND CONCLUSIONS 

With narrow-band optical imaging of seven clusters we have shown 
that Ha and [Nil] emission can be traced to beyond 20 kpc from 
the centre in four clusters previously classified as cooling flows. Al- 
though our images of these sources are deeper than those already in 
the literature, they show very similar morphologies to the previous 
data. A typical non-cooling flow radio galaxy at a similar redshift 
to these clusters shows only nuclear Ha emission when observed 
to the same surface brightness limit. 

There are clear variations in Ha morphology from source to 
source. Two sources are filamentary (Abell 1795, SeiT59-03) and 
two are roughly circularly symmetric (Abell 2204, Abell 2597). 
Abell 2597 shows an azimuthally averaged surface brightness pro- 
file which follows an R~ 3 distribution to beyond 40 kpc (limited 
by the signal-to-noise in our data). The profiles of the other three 
clusters drop more rapidly with radius than this. In particular the 
filaments in Abell 1795 and Serl59-03 appear particularly strongly 
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confined (hence the greater depth of our images show little differ- 
ence from previously published images). We are therefore seeing 
a physical limit on the size of the nebulosity, rather than a limit 
imposed by the signal-to-noise of our images. 

Our spectra show that for each source the H2 emission follows 
the Pact emission closely in spatial distribution, surface brightness 
and dynamics across the nebulosity. Clearly all phases of the nebu- 
losity and the source of excitation must be very intimately mixed on 
scales much smaller than our spatial and spectral resolution. Given 
the lack of dramatic variation in molecular to atomic emission as 
a function of radius, the ratio of physical parameters leading to the 
line ratios vary little over scales of 2 to 20 kpc from the centre of 
the cluster. The consequence of this is that we expect to see the 
molecular gas follow the atomic gas out to well beyond 40 kpc in 
systems such as Abell 2597. This has a minimal consequence on 
the total K— band H2 luminosity because of the rapid dropoff of 
surface brightness with radius as noted above. 

Our detailed and sensitive K— band lo ng slit spectrosc opy of 
three of these sources confirms the results o f Jaffe et alJi2001I> that 
the Pa to H2 line ratios differ radically from those found in star- 
burst galaxies. They are indicative of excitation and heating by a 
spectrum harder than O-stars (or more particularly a "bluer" EUV- 
FUV colour). This data does not allow us to determine whether this 
is caused by hotter stars, a substantial soft tail to the X-ray spectrum 
of the ICM, or some other means. The hardness of the spectrum can 
be independent of the ionization parameter in the gas. 

The difficulty of keeping molecular gas warm by any plausible 
radiative excitation suggests that in this case the heating sources 
must be very inhomogeneous distributed, and in close association 
with the gas clouds. This inhomogeneity of the gas is supported by 
a comparison of CO emission and X-ray absorption measurements 
in at least 2 clusters. 

The velocity dispersion curves of the molecular and ionized 
gas are sharply peaked within 8 kpc of the centres of Abell 2204 
and 2597 with dispersions typically peaking at a ~ 200 — 300 km 
s _1 . All previously published if— band spectra of cooling flows 
are of emission from this region. The motions in this region are 
probably dominated by the energetics of the AGN. In the case of 
SeiT59-03, our A"— band slit traced the large-scale emission fila- 
ment and did not cross the nucleus of the galaxy and avoids this 
high velocity dispersion region. Beyond 8 kpc the velocity disper- 
sion drops below the instrumental resolution of approximately 50 
km s _1 . 

The mean velocity of the gas relative to the centres of each sys- 
tem show values of typically ±100 km s _1 and in Abell 2597 and 
SeiT59-03 show evidence for coherent shear or rotation over 20- 
30 kpc. These velocities are much smaller than the typical stellar 
velocity dispersions in brightest cluster galaxies and are therefore 
insufficient to support the molecular gas against the gravitational 
field of the central galaxy. Moreover, the density of the molecular 
gas indicates that it cannot be pressure supported in equilibrium 
with the hot ICM. Either the molecular gas is continually replen- 
ished from a pressure-supported reservoir of gas on timescales of 
10 Myr or is required to be supported by some other mechanism. 
The kinematics of the molecular gas at large radii seems to exclude 
an origin for the gas either in material stripped from infalling galax- 
ies or material pushed out of the centre of the system by the action 
of the central radio source. 
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